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Abstract 

By analyzing the external field dependence of correlation functions, the magnetic properties of 



hadrons can be determined using lattice QCD in magnetic fields. To compute the magnetic mo- 
ments and polarizabilities of charged hadrons, for example, one requires sufficiently weak magnetic 
fields. Such field strengths, however, lead to closely spaced Landau levels that are not straight- 
forwardly resolved using standard lattice spectroscopy. Focusing on charged spinless hadrons, we 
\ introduce a simple projection technique that can be used to isolate the lowest Landau level. As the 

technique requires the explicit coordinate-space wave-function, we investigate the extent to which 
the continuum, infinite volume wave-function can be employed. We find that, in practice, the ef- 
fects of discretization can be handled using a perturbative expansion about the continuum. Finite 
volume corrections are taken into account by using the discrete magnetic translational invariance 
of the torus. We show that quantized magnetic fields can lead to pernicious volume effects which 



depend on the magnetic flux quantum, rather than on the lattice volume. 
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I. INTRODUCTION 



Quarks are unique in the Standard Model because they are the only particles charged 
under all three gauge groups. While electromagnetic interactions of quarks are dwarfed 
by their strong interactions, electromagnetic observables and tiny electrodynamic effects 
ultimately give one an intuitive picture of the underlying QCD dynamics present in the 
vacuum and in hadrons. Studying QCD in the presence of classical electromagnetic fields, 
moreover, allows one to address how QCD responds to external conditions. In recent years, 
considerable steps have been taken to investigate QCD in external electromagnetic fields 
using lattice gauge theory techniques; for an overview, see [H, 0]. 

The external field method in lattice QCD appeared quite early on. Background magnetic 
fields were used in the first calculations of the nucleon magnetic moments 0, 0]. Soon 
thereafter, it was realized that effects at second order in background electric fields could 
be used to extract the electric polarizabilities of neutral hadrons [5J. In the context of 
weak external magnetic fields, the background field method was revisited several years ago 
to compute magnetic moments and polarizabilities of hadrons 0, 0]. These studies were 
completely quenched, and non-uniform magnetic fields were employed with effects resulting 
from field gradients mitigated by the imposition of spatial Dirichlet boundary conditions. 
Such boundary conditions produce non-perturbative effects that are unfortunately difficult 
to quantify. 

Uniform magnetic fields can be achieved on the lattice by imposing a quantization con- 



dition on magnetic field strength |8Hl0j. In QCD, this condition takes the form 



i in 27m * n\ 
\q d \B = —jy, (1) 

where L is the spatial size of the lattice, which we assume to be the same in each spatial 
direction, and \q&\ = ~e is the magnitude of the electric charge of the down quark. The 
integer n$ is the magnetic flux quantum of the torus. A study of magnetic moments of spin- 
3/2 resonances employing uniform magnetic fields appeared a few years ago [11]. While the 
quantization condition has generally restricted one to prohibitively large magnetic fields, 
lattice volumes are ever increasing, and this trend will consequently lead to lattices that are 
able to support several quantized values of perturbatively small magnetic fields. 

Due to the existence of Landau levels, the correlation functions of charged hadrons de- 
pend strongly on the magnetic field. For a hadron of mass M, the characteristics of the 
Landau levels depend on the dimensionless ratio \eB\/M 2 . When the magnetic field is 
large, \eB\/M 2 3> 1, the Landau levels are widely separated in energy. In this regime, 
one can expect that standard lattice spectroscopy will filter out the lowest Landau level in 
the long Euclidean time limit. On the other hand, for small values of the magnetic field, 
\eB\/M 2 <C 1, the energies of Landau levels are very closely spaced and standard spec- 
troscopic techniques will be of limited use. To deduce the magnetic properties of charged 
hadrons J3 one is ultimately interested in the regime \eB\/M 2 <C 1, however, in the near term 
one may have to settle for more intermediate values, for which \eB\/M 2 < 1. Nonetheless, 
as the available magnetic fields become smaller, the need to explicitly treat Landau levels 
becomes greater. Fortunately this can be achieved easily using the proper-time formalism, 



1 As these properties are dynamically determined, one must additionally worry about energy scales set 
by the virtual degrees of freedom within hadrons. Because the charged pion is the lightest electrically 
charged state in QCD, the ratio \eB\/m% is important for all hadrons. 
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which was pioneered in quantum field theory by Schwinger [12J. In this work, we develop 
the technique for lattice correlation functions by focusing on charged scalar particles. 

In an external magnetic field with \eB\/M 2 < 1, the tower of Landau levels complicates 
standard lattice spectroscopy. As we show, one need not rely upon the long-time limit 
of the correlation function to filter out the lowest Landau level. One can project out this 
state explicitly, and thereby considerably reduce the challenges inherent to studying charged 
hadrons in external magnetic fields. We begin in Sec. [TT] by detailing the method using a 
continuum action defined in infinite volume. For simplicity, we focus on charged scalar 
particles. To address the effects of discretization, we consider a charged particle effective 
action on an infinite lattice in Sec. IIHI We find that the effects of discretization can be treated 
in practice using perturbation theory with a continuum wave-function. Next in Sec. IIVt we 
investigate finite volume corrections. Using the discrete magnetic translational invariance 
of the torus, we account for finite volume corrections to the lowest Landau level. While 
volume corrections are generally suppressed by exponentially small terms, we expose certain 
corrections that are not exponentially small, and discuss how the infinite volume limit can 
be recovered. This subtlety owes to the magnetic field quantization required on the torus. 
Finally we conclude in Sec. |V]with a brief summary of our findings, and directions for future 
work. 



II. CONTINUUM METHOD 

The object for studying the spectrum of QCD is the two-point correlation function. Stan- 
dard lattice spectroscopy relies on the large Euclidean time limit of the two-point correlator. 
With a mixed momentum-time representation, the two-point correlator is^| 

G(p,r) = J dxe ip - s (0\®(x, t)&(0,0)\0), (2) 

where for simplicity we restrict our attention to a charged scalar hadron which has an 
interpolating field <3>. The Fourier transform introduces the coordinate-space wavefunction, 
ipp{x) = e ip ' x , of the momentum eigenstate p. The spectral decomposition of the correlator 
shows that the ground state saturates the correlation function in the long Euclidean time 
limit, with excited-state contamination suppressed by the factor ~ exp(— AEt). Here AE 
is the splitting between the ground and first excited states, and we have assumed that the 
overlap factors for both states are approximately equal. For Euclidean times on the order of 
r ~ 1 fm, a natural energy-level splitting, AE ~ Aq CD , ensures reasonably good suppression 
of the excited states. 

The addition of a magnetic field, B, considerably complicates the spectral decomposition 
of the two-point function. Energy levels of multiparticle states, for example, become cum- 
bersome to enumerate. We will assume the zero-field ground state is the single-hadron state 
<p. This ground-state hadron, which is a state possessing the quantum numbers of $, now 
consists of an infinite tower of Landau levels. These energy levels are not widely separated in 
weak magnetic fields. For an ideally weak external field, the energy splitting between neigh- 
boring Landau levels is set by AE = \eB\/M, where M is the hadron's mass. In a given 

2 Throughout we omit explicitly writing the end-points of integration for all integrals spanning the entire 
real line. Thus we write J dx for the integral f_°° dx. Similarly for sums over all integers, we do not 
explicitly write out the limits. Thus we use the notation ^ n for the sum Yln=-oo- 
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magnetic field, the proton's Landau levels, for example, will be more closely spaced than the 
charged pion's leading to comparatively more challenging lattice spectroscopy. Even for the 
best case, the pion itself, the energy level splitting in perturbatively weak fields requires r 
considerably greater than 1 f m to resolve. Consequently it is beneficial to handle the Landau 
levels explicitly. 

The single-particle effective action for the composite field in a uniform electromagnetic 
field has the form 

S = J d 4 x [D^D^ + MV + c F^F^(f> + c 2 F p{ ^F u}p D^D u (fl , (3) 

which follows on account of gauge invariance, parity invariance, and Euclidean invariance. 
The gauge covariant derivative is defined in the usual way, D^cf) = d^cp + ieA^cj), and the 
electromagnetic field-strength tensor is F^ v = d^A u — d u A^. The curly braces denote the 
symmetric traceless part of a tensor, namely 0{ pu } = | (O pu + O vpi — ^5^ v O pp ). In the 
effective action, we have written down all terms containing up to two powers of the field- 
strength tensor. There are of course higher-order terms allowed by symmetries, but we 
assume the external field to be sufficiently weak. The low-energy constants, c and c 2 , 
describe the helicity, A = and A = 2, couplings of two photons to the (p. From a simple 
matching calculation, these low-energy constants can be identified as 

Atx 

c = 7iM(a E - M ), c 2 = - — (a E + f3 M ), (4) 

where a E is the electric polarizability of the 0, and /3m its magnetic polarizability. 
Now we specialize to the case of a magnetic field specified by the vector potential! 

A, = (-5x 2 ,0,0,0). ( 5 ) 
After a field redefinition, the effective action can be written in the form 

S = J d A x [D^D^ + M 2 (t> ] (f\ , (6) 

where M. is given by 

M{B) = M- Unf3 M B 2 + C(-B 4 ), (7) 

and the higher-order terms in the magnetic field that were neglected above are now effectively 
subsumed into the 5-dependence of M.. The goal of a background field lattice computation 
is to determine the magnetic polarizability (3m from the magnetic field dependence of Ai. 
In an external magnetic field, one can consider the two-point correlator, 

G b (t) = J dx (0|$(f,r)$t(o,0)|0) B , (8) 

where the subscripts denote the external field dependence. Because the 22-component of 
momentum is not a good quantum number in the external magnetic field, the Fourier trans- 
form implicit in Eq. ([8]) receives contributions from an infinte tower of Landau levels. Using 



3 Ordinarily one might average over all gauge equivalent vector potentials, however, this is not what one 
does when adding the external field in a lattice QCD calculation. In this situation, the gauge for the 
external field is fixed, and charged particle correlators are gauge variant. 



4 



the effective action in Eq. (J6]) and following 13j, the contribution to the correlator can be 
written in the forrrfl 

G b (t) = / da (9) 
Jo a/ s cosh [ens ) 

where the non-standard r-dependence arises from summing contributions from the entire 
tower of Landau levels. In the extreme long-time and non-relativistic limits, only the lowest 
Landau level contributes, Gb{t) ~ e-( M +2\ eB \) T , As lattice calculations are far from these 
ideal limits, one must confront contributions from higher Landau levels. The explicit form 
of the correlation function in Eq. (Q could be used, in principle, to fit the lattice correlator^ 
however, there are hidden assumptions. In deriving this expression for the correlator, we 
integrated over all space. This leads to contributions from Landau levels with an average 
size extending a considerable distance from the origin. Such contributions are modified due 
to boundary conditions. It is thus highly desirable to filter out only the lowest Landau level 
which has the smallest average size. 

To remedy the situation with narrowly separated energy levels, we return to the correla- 
tion function in Eq. (jSJ). By using this form, one is implicitly assuming the single-particle 
wave-function, ip* (a^) = e tp2X2 , of a free particle with p 2 = OU cf. Eq. ([2]). A more desirable 



Gb(t) = / d$i$°\x 2 ) (0|$(£,r)$t(0,0)|0>B, (10) 



P2 

correlation function is 



where i/)q(x) is the ground-state harmonic oscillator wave-function 

4°)( x ) = e-^ 2 . (11) 

The absolute normalization of this wave-function is irrelevant for our purposes, and the 
notation for the wave-function is chosen to be consistent with Sec. II III below. Computing 
the 4> contribution to the correlator in Eq. ({TO]) produces a simple result. Using Schwinger's 
proper-time trick, we find an exponential fall-off of the correlator 

g B (r) = Z^e- EoT , (12) 

with the energy given by 

E = VM 2 +\eB\. (13) 

The new correlation function Gb{t) allows one clean access to Ai, and consequently to the 
magnetic polarizability (3m- 



4 The overlap factor, which is generically written as Z^,, also contains kinematic factors. We will write 
the same overlap factor throughout; although, it may have different numerical values in different 
correlators. Because the overlap of the interpolating field 4> with the single particle state is a priori 
unknown, we do not keep track of such kinematic factors. 

5 The utility of external field correlators in fitting lattice data was originally suggested for background 



electric field computations 14]. Analogous to the Landau level problem, charged particle two-point 



functions in external electric fields have non-standard r-dependence, and the predicted functions have 



been sucessfully used to extract properties of pseudo-scalar mesons and baryons [la, Il6|. 
6 If one imposes spatial Dirichlet boundary conditions, there is an analogous issue. The sum over spatial 
lattice sites produces a tower of standing waves. This tower can be directly avoided by projecting out the 
lowest mode using its coordinate wavefunction, vp(x) = sin(^). While this method would circumvent 
the momentum problem, it unfortunately does nothing to resolve the myriad difficulties with Dirichlet 
boundary conditions, e.g. U(\)a restoration, the possibility of chiral symmetry restoration, etc. 
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An obvious drawback of the projection method is that the coordinate-space wave- function 
for the lowest Landau level is required. On a finite space-time lattice, this wave-function 
does not have the simple GauBian form well-known from elementary quantum mechanics. 
Modifications to the wave-function will arise from the discretization of space, as well as the 
finite spatial volume. We consider each of these effects below. 



III. DISCRETIZATION EFFECTS 



We investigate first the effects of discretization on the lowest Landau level. To this end, we 
temporarily ignore the boundary conditions. Addressing finite-size effects is postponed until 
the following section. In the absence of boundary conditions, a natural question remains. 
Can one reasonably approximate the lowest Landau level using a continuum wave-function? 
To investigate the answer to this question, we consider the theory on an infinite spatial 
lattice with a uniform lattice spacing a. The lattice sites are labeled by a vector of integers, 
n = (ni, n 2 ,n 3 ), which corresponds to the coordinate vector x = an. For simplicity, we keep 
the time direction continuous. 

The projection method obviously generalizes to the lattice, where one now requires the 
lowest lattice Landau level. It is described by a wave-function, where we use the 

subscript n 2 to index the lattice-site dependence. Generally we use subscripts for discrete 
indices, for example, the lattice version of the interpolating field is now written $^(£4). 
With the wave-function, if)o tna , we then form the discrete analogue of Eq. (JT0|) . 



^(r) = ^^ in2 (0|$,(r)$t(0)|0) 



B- 



(14) 



Consequently the sum over n 2 in the lattice correlator, Eq. (I14p . will project out only the 
lowest energy state in the tower of lattice Landau levels. 

To determine the wave-function for the lowest lattice Landau level, we must know the 
lattice form of the single-particle effective action, which is the discretized version of Eq. 
Using this action, one solves the discrete eigenvalue equation for the lowest energy eigenstate, 
and ipo t n 2 is then the corresponding eigenvector. We will assume a highly plausible form for 
the discrete action of the charged scalar. When one couples gauge fields to matter fields 
on a lattice, gauge invariance requires the use of link variables, U^ x = exp(ieaA^ x ). A 
lattice form of the effective action consistent with cubic symmetry, gauge invariance, and 
Hermiticity appears as 



5 = a 3 ^ J dx4, (fyix^Vfttf <t>n'{x4), 



with 



d 2 



Af 



(15) 



(16) 



In writing Eq. ( Tl6|) . we have made use of the vanishing of the temporal component of the 
gauge field, and the time-independence of the gauge field. In the absence of the external 
field, Eq. f ll6p leads to a dispersion relation of the form 

3 



E 2 = M 2 + 4S sin2 ( a ^/ 2 )' 



(17) 
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One can use the measured spatial momentum dependence of scalar particle energies to 
establish the precise form of the lattice dispersion relation in a given lattice calculation. For 
our present purposes, we assume the dispersion relation to be given by Eq. ([IT]) , and the 
discrete action in Eq. (JTBjl is merely the corresponding gauged action. 

Returning to the assumed form for the discretized single-particle action in Eq. f fl6|) . we 
see that the gauge links for the external field, Uj^, introduce explicit lattice-site dependence. 
As a result, the ^-component of lattice momentum is not a good quantum number. This 
mirrors the situation in the continuum. The remaining spatial lattice momenta, pi and p%, 
are good quantum numbers, as is the time-component of momentum, p 4 . It is thus useful 
to take the Fourier transform; namely, in terms of p^ = (pi, 0,Ps,P4), we have 



4>nM = Z) / dxte^M^)- (18) 

ni,ri3 

In terms of the Fourier transformed field, the action now appears as 

ri2,n' 2 a 

Because we can consider correlation functions with two vanishing components of spatial 
lattice momentum, p± = p% = 0, we restrict our attention to this sector, in which 

V n2 , n > 2 (p A ) = T n2 y 2 + 5 n2 y 2 [pi + M 2 + V(n 2 )} , (20) 

where the discretized kinetic term is 

T n2j n' 2 — J [3ri2,ri' 2 +l + $n 2 +l,n 2 ~ 25 n2>n ^] , (21) 

and the potential term is 

V(n 2 ) = 4 sin 2 (ea 2 Bn 2 /2) . (22) 
a 2 

The eigenfunctions that diagonalize the operator T + V are solutions to the lattice Landau 
level problem. The lowest level, for example, satisfies the coordinate-space eigenvalue equa- 
tion, (T + V)ipo^ n2 = Ao^o,n 2 - As a consequence, the contribution to the lattice two-point 
correlation function in Eq. fll4p has the simple exponential form 

Q b {t) = V~ £ ° T , (23) 

with the energy eigenvalue given by 



So = VM 2 + V (24) 

Let us investigate the form of the eigenvalue equation for the lattice Landau level problem 
more closely. In the auxiliary proper-time quantum mechanics, the operator T + V appears 
as (twice) the Hamiltonian, and can be written in the abstract form 



4 

T+V = — 

a 2 



sin 2 (ap 2 /2) + sin 2 (eaBx 2 /2) . (25) 
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In the continuum limit, we have (twice) the Hamiltonian of the simple harmonic oscillator 

T + V =p 2 2 + e 2 B 2 x 2 2 , (26) 

which respects the canonical transformation (£2,^2) |eS|x2^. The lattice Landau 

level problem, Eq. ( |25|) . also respects this canonical transformation. 

The dimensionless parameter governing the expansion of the Hamiltonian in Eq. ( 125|) 
about the continuum limit is the magnetic field in lattice units, namely 

b=\ea 2 B\. (27) 

As quantized magnetic fields on a lattice will satisfy Eq. (TI]), we have b = 6ir\n<s>\/N 2 , where 
N s is the number of spatial sites. Typically one needs at least four values of the magnetic 
field strength to extract a polarizability with confidence, so (n$) max = 4. The number 
of spatial sites must be at least N s = 32 to have pertubatively small fields compared to 
hadronic scales with a typical lattice spacing of a = 0.1 fm. Thus in practice, we are limited 
to b < 1/10, and an expansion in b is well justified. 

When we are near the continuum limit, the wave-function of the lowest Landau level, 
, 4 , o > \x2) in Eq. ffTTT) . restricts the coordinate to values X2 < |e.B|~i. Thus when we series 
expand the potential V about the continuum limit, V = YljLo ^2?' ^ ne ^"^ n t erm 

V 2j oc\(eaBx 2 ) 2 i +2 <VV , (28) 
cr 

is suppressed by j powers of b. The kinetic operator T can also be expanded near the 
continuum limit, T = YlJLo^j- In the lowest Landau level, the momentum p 2 is also 

restricted, P2 < \eB\^ on account of the canonical transformation. Consequently, the j-th 
term in the expansion 

T2^- 2 {ap 2 ) 2 ^ 2 <VT Q , (29) 

is also suppressed by j powers of b. Hence we address the leading-order effects of discretiza- 
tion by treating T 2 + V 2 as a perturbation of the continuum result. 
The explicit form of the expansion to second order is 

T 2 + V 2 = -j^ [(ap 2 ) A + 6 4 {x 2 /af] . (30) 

The computation now amounts to basic Rayleigh-Schrodinger perturbation theory. We com- 
pute the leading correction to the energy eigenvalue, as well as the wave-function. Writing 
the eigenvalue of the lowest lattice Landau level as 



A = Ai 0) + A( 1) + ---, (31) 



with a 2 Ag°' ) = 6, we have at once 



a 2 A W = _° (32 ) 



Having deduced this correction, we can write down the perturbative expansion in b of the 
the energy eigenvalue Sq in Eq. f )24p . Given in lattice units, we find 

a 2 S 2 = a 2 M 2 + b - Q + 0) b 2 + 0(b 3 ), (33) 
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FIG. 1. Discretization corrections to the coordinate wave-function of the lowest Landau level. 
Shown as a function of ri2 is the relative difference Aipo(ari2) defined in Eq. ([57|) . In lattice unites, 
the magnetic field strength is taken to be b = 0.1. Values this size and smaller are required in 
practice. 

where the /3 term depends on the magnetic polarizability 0m in the form 

13 = aM^, (34) 
a a 6 

with a = t~. Polarizabilities of hadrons have values of natural size in units of 10 _4 fm 3 , 

4-7T ' 

although one expects smaller values for polarizabilities at larger-than-physical pion masses. 
For a physical mass pion, M = m n , and a lattice spacing of a = 0.1 fm, we expect /3 ~ 1. 
From this value, we see that the discretization correction to the energy of the lowest Landau 
level will affect the extracted value of the polarizability by ~ 10-20% if ignoredQ 

To investigate the discretization corrections to the ground-state wave-function, we simi- 
larly write the coordinate wave-function as a perturbative expansion about the continuum 
limit 

M^)=^o\x2)+4 1) ^2) + -- - , (35) 

and we find 

^\x 2 ) = ^=^\x 2 ). (36) 

For magnetic field strengths of size b < 0.1, the correction to the wave-function ip^\x2) in 
Eq. (1361) is quite small. We can exhibit this by considering the difference in the ground-state 
wave-function relative to the continuum one, 

AMx2) = ^y— • (37) 

% '{x 2 ) 

This relative difference is plotted in Fig. [T]for the value b = 0.1. In practice, magnetic fields 
of this size and smaller are required to probe the magnetic polarizability of the <fi. From the 
figure, we see that corrections to the continuum wave-function will be less than 1% for such 
magnetic fields. 



7 With respect to finite lattice spacing effects, we must also remark that the electromagnetic current 
receives 0(a 2 ) corrections. Such discretization corrections have nothing to do with Landau levels, and are 
thus present even for neutral hadron correlation functions. Discretization corrections to the polarizability 
due to the electromagnetic current are proportional to (oAqcd) 2 ; which, on a lattice with a = 0.1 fm, is 
expected to be < 4% for Aq C d < 400 MeV. 



9 



Notice that the relative correction actually grows as a function of n 2 because it is propor- 
tional to a Hermite polynomial, namely H^y/bn^). While this polynomial is unbounded, 
large corrections only occur when the wave-function is exponentially small. For example, 
to get a 10% discretization correction to the wave-function in a field of b = 0.1, one must 
go out to a distance n 2 ~ 9 from the origin. But the ratio of the wave-function at this 
point to its value at the origin is ip < ^\an 2 ) / ipf\^) ~ 10 -2 . Because a determination of the 
correlation function itself to much better than 1% is not likely, the resulting discretization 
correction to the wave-function can be safely neglected. Said another way, plots of iPq\x 2 ) 
and ip^\x 2 ) + i/}q(x2) for b = 0.1 are indistinguishable, unless we zoom in on the tails of 
the GauBian. 

In practice, one will need to verify the lattice form of the dispersion relation for the particle 
of interest in order to address the effects of discretization. While the form of the dispersion 
relation need not be exactly that assumed in Eq. (1T71) . our power-counting arguments readily 
generalize to other forms. For small values of b, only the 0(a 2 ) corrections to the proper-time 
Hamiltonian will be required 

AH = ~lL? \- Cl ^ p2)A + C ^\x 2 /af] . (38) 

This differs from Eq. (|30|) by the introduction of two arbitrary coefficients, C\ and C 2 . The 
coefficient C\ can be determined numerically by studying the deviation of the dispersion 
relation from the continuum one, 

£ 2 = M 2 +p 2 -^X>) 4 + "--, ( 39 ) 

where we have made use of cubic invariance, and eliminated the term proportional to (p 2 ) 2 
by making a field redefinition. The coefficient C 2 appearing in Eq. fl3"8l is not an independent 
parameter; it is fixed by gauge invariance. The effective action describing a particle with 
Eq. f[3T?j) as its dispersion relation is gauged by the replacement pj — > pj + eAj. After 
projecting the good components of momentum, p\ and p^, to zero in such an effective action, 
we arrive at Eq. (|38|) with C 2 = C\. Notice that additional terms, such as (p 2 ) 2 (x 2 ) 2 , are 
forbidden by the combination of cubic symmetry and gauge invariance. Thus to address the 
effects of discretization on the lowest Landau level, we perturb about the continuum with 
the interaction 

AH = -^- 2 [(ap 2 r + b\x 2 /ar}. (40) 

The results presented above employ the value C\ — 1, however, this coefficient should be 
determined non-perturbatively for the particular lattice QCD action employed in actual 
calculations. For the pseudoscalar Goldstone mesons, one generally expects the breaking 
of rotational invariance in Eq. ( 13 9 p to be quite small (l7| . while for baryons this is not 
necessarily the case [l~8|. 



IV. FINITE VOLUME EFFECTS 

Having established that discretization corrections can be treated in a perturbative ex- 
pansion about the continuum, we now turn to finite volume effects. To study these effects, 
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we consider the theory defined on a continuous torus of length L in each of the three spatial 
directions, with the temporal extent kept infinite. We need not consider a discrete torus, 
i.e. a finite lattice, because the combination of both discretization and volume corrections 
is a doubly small effect. To determine the finite volume modifications to the projection 
technique, we must understand the behavior of charged particle two-point functions under 
spatial translations by L. 

On a torus, there is a discrete magnetic translation group, see e.g. [HI], a f ew aspects of 
which are required to address volume corrections. The gauge potential given in Eq. (jSJ) 
is periodic up to a gauge transformation 

Ax{x + x 2 L) = A x {x) + <9iA(x), (41) 

with all other directions strictly periodic. The gauge transformation function appearing 
above is A(x) = —BLx\. The corresponding gauge transformed matter field must then 
satisfy the magnetic periodic boundary condition 

<f){x + x 2 L) = e ieBLxi <f){x), (42) 

with the other directions periodic. The boundary conditions in the x% and x 2 directions 
can only be consistent if the magnetic flux through the X\, a^-plane is quantized. We will 
label coordinates in this plane simply by x±. For down quark fields of fractional electron 
charge, the quantization condition is given in Eq. (pQ); and, in turn, this gives a magnetic 
flux quantum of N& = 3n$ for a scalar hadron having unit charge. 

To generalize the projection technique to the torus, we must isolate the lowest Landau 
level. The ground-state wave-function on a torus must be a sum of images of the infinite 
volume wave-function [l9| 

VT(^) = E + uL)e-^ iN *^ L , (43) 

V 

in order to maintain magnetic periodicity. Notice that this particular finite volume wave- 
function corresponds to a state with zero momentum in the Xi-direction. With a lattice 
determination of the two-point functiorJl 

G f b v (x ± ,t) = y^a:3(0|$(f ± ,X3,r)<l>t(0 ± ,0,0)|0)^ (44) 

we then form the finite volume generalization of Eq. ffTO]) 

Gb V (t) = J l dx ± r FV (x ± ) G™(x ± , T ). (45) 

Appealing to Schwinger's proper-time trick, and a resolution of the identity in the infinte 
xj_-plane in terms of cells of area L 2 indexed by integers {yi, u 2 ), we find the <fi contribution 
to the correlator indeed has a simple form 

Q F B V {r) = Z^e- E -\ (46) 

8 We consider the lattice to extend from — U to in each of the spatial directions. This is an inessential 
choice, and serves to simplify the discussion of finite volume corrections in the winding number expansion. 
For a lattice that extends from to L in each of the spatial directions, the relative change in the ground- 
state wave- function is exactly the same as in Eq. (T15)) . As X2 now extends to L in this case, there is a 
very large finite volume correction from the v = 1 image. This is essentially the peak of the ground-state 
wave-function located at the origin, because, due to magnetic periodicity, the origin, X2 = 0, is equivalent 
to X2 = L up to an Xi-dependent phase. 



11 



«2 



FIG. 2. Finite volume corrections to the coordinate wave-function of the lowest Landau level. 
The plot shows the relative difference Ai/jq V (xi = 0,^2 = an2) defined in Eq. (|48|) as a function 
of n%. We choose the minimal flux quantum, N$ = 3, and a relatively small lattice size, L = 24a, 
for which the corresponding value of the magnetic field is b = 0.033. Finite volume corrections for 
x\ ^ are bounded in absolute value by those shown in the plot. 

with the energy, Eq, the same as in infinite volume, Eq. (|T3l) FI 

Taking a pion interpolating field for $, the coordinate-space correlator in Eq. (144 p remains 
positive definite in an external magnetic field. In projecting out the ground state with the 
wave- function in Eq. P3|) . we thus restrict our attention to the real part of Q^ v (r), which 
depends on the real part of ipQ V {x±). Consider the change in the finite volume wave-function 
relative to the infinite volume one, 

A ^ (fJ = Wfc) . (47) 

This relative change takes the form 

oo 

A^ v (x ± ) = 2^2 e~ M * cos(2ir NqUXx/L) cosh(2vr^z/x 2 /L). (48) 

v=\ 

Physically we expect volume corrections to be largest in the smallest magnetic field. The 
smallest magnetic field allowed on the torus has the flux quantum = 3. A plot of the 
relative change in the wave-function AipQ (x±) is shown in Fig. [2] for a magnetic field with 
flux quantum iV$ = 3 on a lattice of size L = 24a. The relative change shown in the figure 
overwhelmingly arises from the winding number v = 1. Contributions from v = 2, for 
example, are suppressed by at least e _67r , as we show below. Not surprisingly, the largest 
finite volume corrections arise at the ^-boundary of the lattice. 

Let us carefully analyze the relative change Ai\)§ v {x±) by considering the contribution 
from a generic winding number v in the series. The Xi-dependence of this contribution enters 
through the cosine, cos (2iiN^i'Xi/L), and so is always bounded by unity. The z/ th term in the 
series depends exponentially on the magnitude of x<i through the hyperbolic-cosine. Because 
the x 2 -coordinate is bounded, |x 2 | < ■§, however, we have cosh (27iN< s> b / X2/ L) < |e y7riV * to 
exponential accuracy. Thus the contribution to the finite volume effect from a generic 



9 It should be noted that there are also dynamical finite volume corrections which affect M 2 that we are 
neglecting. These corrections are suppressed by exp(— m„L), however, their treatment is subtle due to 



the holonomy of the external field 2CJ, [21 1 
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winding number v is bounded by e - u ( u ~ 1 ) 7rN <s> ; i. e . finite volume effects are suppressed by an 
exponential factor involving the magnetic flux quantum. This suppression occurs provided 
the winding number is greater than one. When v = 1, there is no exponential suppression 
in the relative change when x 2 approaches the boundary of the lattice. 
At the ^-boundary of the lattice, we have 

A^ v ( Xl = 0,x 2 = ±L/2) = 1, (49) 

to exponential accuracy, which is independent of the flux quantum N$. While this corre- 
sponds to the maximal value for the relative change, we should also consider the size of the 
ground-state wave-function at the ^-boundary, because this is where the wave-function is 
the smallest. The ratio of the infinite volume wave-function at the boundary to its value at 
the origin is given by 

4 0) (x 2 = ±L/2)/^ 0) (0) = e-f^, (50) 

which is remarkably independent of the lattice size L. Unlike discretization corrections to 
the wave-function, finite volume corrections are generally non-negligible. For example, in the 
smallest magnetic field, N<& = 3, the boundary-to-origin ratio ~ 0.1 is not considerably less 
than unity and one must consider the ~ 100% correction to the infinite volume wave-function 
from the v — 1 image. 

The independence of this result on the lattice size L brings us to our final point. There is a 
subtlety in taking the infinite volume limit. The finite volume correction in Eq. (I48p does not 
vanish in the limit L — > 00. This is due to the magnetic field quantization condition. As L is 
made larger, the magnetic fields allowed on the torus becomes smaller. When the magnetic 
field becomes smaller, however, the ground-state wave-function spreads out in space. In 
this way, increasing the volume need not result in smaller finite volume effects. The tacit 
assumption in this approach to infinite volume is that the flux quantum remains fixed 
as L —7- 00. This is likely to be the case in practice, as one will generate larger lattices and 
then focus on the smallest allowed values of the magnetic field on these lattices. To recover 
the infinite volume limit, however, the magnetic field should be fixecf^l as L — > 00, which in 
turn requires iV$ cx L 2 . While there are still ~ 100% finite volume corrections to the wave- 
function from Eq. (jUJ), these occur at the ^-boundary which is now precisely where the 
wave-function is exponentially suppressed in the volume, cx e~s^ eB ^ L . With such corrections 
buried in the tails of the GauBian wave-function, the infinite volume limit is recovered. 

V. SUMMARY 

Above we explore the correlation functions of charged spinless hadrons in external mag- 
netic fields. For magnetic fields that are small compared to the hadron's mass, \eB\/M 2 <C 1, 
the hadron's magnetic properties can, in principle, be measured from lattice QCD simula- 
tions. Such computations, in practice, require modification of standard lattice spectroscopic 
techniques. This modification is necessitated by closely spaced Landau levels that cannot 
be cleanly resolved from the long-time limit of Euclidean correlation functions. To handle 

10 Actually one need not be this restrictive to recover the infinite volume limit. As long as the flux quantum 
scales with some positive power p of the lattice size, the infinite volume limit will be recovered. The 
finite volume correction to the wave- function remains ~ 100% at the a^-boundary, but the size of the 
wave- function at the X2 -boundary is oc exp (— ■jN®) which will be exponentially small in the volume for 
N$ cx LP. 
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this complication, we develop a projection technique to isolate the lowest Landau level. The 
technique requires a modified two-point correlation function, Eq. (fTUj) . that depends on the 
coordinate-space wave-function of the lowest Landau level. To put the technique into prac- 
tice, one needs to write out the coordinate-time dependence of the correlator data. This 
procedure allows one to later convolve the lattice data with the ground-state wave-function, 
as one may wish to use a continuum wave-function, or include effects from discretization 
and the finite volume. 

We investigate the effects of discretization on the lowest Landau level by assuming a 
form for the discrete action of charged scalar field. For typical lattice sizes, the effects 
of discretization can be treated in perturbation theory about the continuum limit. The 
perturbative correction to the eigenvalue and wave-function of the lowest Landau level are 
given in Eqs. f )32|) and fl36|) . respectively. The discretization correction to the energy could 
affect the extraction of the polarizability at the ~ 10-20% level. On the other hand, the 
discretization corrections to the wave-function of the lowest Landau level are shown to 
be negligible. Effects of boundary conditions are also considered. The wave-function for 
the lowest Landau level on a torus, Eq. (f4"3"j) . is used to account for finite volume effects. 
While the size of finite volume corrections is generally set by an exponentially small factor 
involving the magnetic flux quantum and the winding number is, there is no exponential 
suppression from contributions with v — 1. We show that such finite volume corrections 
near the lattice boundary can be important even as the lattice volume increases. This 
pernicious effect owes to the magnetic flux quantization required on a torus. Caution must 
be exercised in assessing the size of finite volume corrections to the wave-function even 
on large lattices. Nonetheless, we show how to determine the discretization and volume 
corrections to the projection technique, and thereby how to isolate the lowest Landau level 
from lattice correlation function data. 

Further work will allow one to extend the technique. Inclusion of spin degrees of freedom is 
necessary to handle proton correlation functions. As the typical magnetic field strengths are 
semi-relativistic, one must account for the fact that the current and anomalous magnetic 
moment operators in the effective proton Hamiltonian do not commute. This puts into 
question the method of using Zeeman splittings to determine magnetic moments. Another 
avenue for investigation is to consider the effects of discretization in external electric fields. 
The functional form of single-particle correlation functions employed by [l5|, 16[ is modified 



by these effects. Fortunately these studies employed anisotropic lattices with a temporal 
lattice spacing of a t = 0.035 fm. On isotropic lattices, however, discretization effects might 
be important. One can use the methods established here to investigate the finite lattice 
spacing corrections to the extraction of the electric properties of hadrons using background 
field correlators. Finally, although this work deals with weak magnetic fields, there has been 
considerable interest in studying lattice QCD in strong magnetic fields due to applications 



relevant for heavy-ion collisions, see, for example, [22H29| . As stronger magnetic fields are 
more susceptible to discretization effects, one could use the methods developed here to 
investigate finite lattice spacing corrections in chiral perturbation theory with strong-field 
power counting, \eB\ ~ m%, or in perturbative QCD. 

Sustained progress in lattice QCD has lead, in particular, to the generation of gauge 
ensembles on large physical volumes. Such volumes will soon permit the study of QCD 
properties in uniform magnetic fields that are perturbatively weak compared to hadronic 
scales. This opens up the possibility to study the magnetic properties of hadrons using the 
external field method. With the projection technique developed here, charged hadrons can 
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be studied by confronting the Landau levels directly. We look forward to exploratory lattice 
studies employing the technique. 
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